Numerical Results for Ground States of Spin Glasses on Bethe Lattices 
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The average ground state energy and entropy for ±J spin glasses on Bethe lattices of connectivi- 
ties k + 1 = 3 . . . , 26 at T = are approximated numerically. To obtain sufficient accuracy for large 
system sizes (up to n = 2048), the Extremal Optimization heuristic is employed which provides high- 
quality results not only for the ground state energies per spin e^+i but also for their entropies Sk+i- 
The results show considerable quantitative differences between lattices of even and odd connectivi- 
ties. The results for the ground state energies compare very well with recent one-step replica sym- 
metry breaking calculations. These energies can be scaled for all even connectivities k + 1 to within 
a fraction of a percent onto a simple functional form, e^+i = EsuVk + 1 — (2Esk + y/2)/y/k + f , 
where Esk = —0.7633 is the ground state energy for the broken replica symmetry in the Sherrington- 
Kirkpatrick model. But this form is in conflict with perturbative calculations at large k + 1 , which 
do not distinguish between even and odd connectivities. We find non-zero entropies st+i at small 
connectivities. While Sk+i seems to vanish asymptotically with l/(k + 1) for even connectivities, it 
is indistinguishable from zero already for odd k + 1 > 9. 
PACS number(s): 75.fO.Nr, 02.60.Pn, 89.75.-k, 05.f0.-a 
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I. INTRODUCTION 

In this paper we study the ground state (T = 0) prop- 
erties of ± J spin glasses on k + I-Bethe lattices M. The 
Bethe lattices in this case are r-regular graphs M with 
r = k + 1. These are randomly connected graphs consist- 
ing of n vertices in which each vertex has a fixed connec- 
tivity of fc + 1. This constraint contrasts with "random 
graphs" g, | in which pairs of vertices are randomly 
connected, leading to a Poissonian distribution of con- 
nectivities around a mean of (c); these graphs will be 
studied numerically elsewhere ||. We explore the large- 
n regime of low-connectivity graphs, k + 1 = 3, . . . , 26, 
which are of great theoretical interest as finite-connected, 
mean-field models for low-dimensional lattice spin glasses 
H E|. A great number of studies have focused on var- 
ious aspects of this conceptually simple model to hone 
the complex mathematical techniques required to treat 
disordered syste ms 111 |7j , H p| , y, [l| or optimization 
problems @ M, |16|, |17], |§ |lg]. In this paper we will try 
to provide a independent numerical check on the validity 
and accuracy of those techniques. 

Our results, in turn, reflect on the flexibility of the 
extremal optimization (EO) heuristic [^0| in find- 
ing approximate but high-quality solutions for ground 
states of spin glasses on an arbitrary graphical struc- 
ture in a reasonable computational time. These are of- 
ten NP-hard optimization problems which are believed 
to require a computational effort that rises faster than 
any power of n to obtain provably exact solutions [p2[ . 
Thus, exact methods as of yet are not able to provide 
results for large- n problems, at least not with significant 
statistics |23| , except for some special cases . 
Furthermore, there is only a small number of capable ap- 
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proximate algorithms available for the study of T = 
properties of spin glasses (2?], |2^, mostly restricted 
to d-dimensional lattice models, and EO provides a dis- 
tinct alternative which will increase the confidence in the 
numerical results available. In previous papers, we have 
demonstrated the capabilities of EO in determining near- 
optimal solutions by reproducing existing results for 3d 
and Ad spin glasses and obtaining new results for the col- 
oring problem pl|, |3C| | and the graph partitioning problem 
|||. The results in this paper show that EO is 
not only capable of approximating ground states well but 
also of sweeping the entire configuration space efficiently 
to determine the degeneracy of ground states |3^] . Un- 
like other methods, EO never "freezes" into local minima 
and proves to be limited only by the ability to store new 
ground states. 

We find that our results for the ground state ener- 
gies are consistent with the theoretical results of the one- 
step replica symmetry-broken (1RSB) solution of ± J spin 
glasses on k + 1-connected Bethe lattices. Our numerical 
result for fc+1 = 3 below clearly excludes the replica sym- 
metric (RS) solution and are consistent with the 1RSB 
results Beyond that our results suggest subtle dif- 
ferences between even and odd values of the connectivity 
fc+1, with no obvious way to continue smoothly between 
them j33| . These oscillations may doom perturbative cal- 
culation for k + 1 — » co H |llj|. While the expectation 
has been raised that the entropy per spin of the ground 
states in this model should be vanishing for any fc + 1 
we find that the entropy is finite and decaying like 
l/(fc + 1) for large, even fc + 1. For odd fc + 1 it is non- 
zero only for small values and may be vanishing already 
beyond some finite, odd connectivity. 

In the following we first introduce the Bethe lattices we 
used in the numerical calculations. In Sec. Ill we briefly 
describe the EO algorithm which is amply discussed else- 
where [^0| ^l], 34j . In Sec. IV we present a few simulations 
to reproduce known results to gauge our procedure. In 
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Sec. we present our nume rical results, followed by an 
extensive discussion in Sec. V C . Some conclusions are 



presented in Sec. VI 



II. SPIN GLASSES ON BETHE LATTICES 

Disordered spin systems on random graphs have been 
investigated as mean-field models for low-dimensional 
spin glasses or optimization problems, since variables are 
long-range connected yet have a small number of neigh- 
bors. Particularly simple are Bethe lattices of connectiv- 
ity k + 1 jy, [j], HI, also called fixed-valence or r-regular 
random graphs |p| [l9| [h], [55|. These are graphs consist- 
ing of n vertices where each vertex possesses a fixed num- 
ber k + 1 of bonds with randomly selected other vertices. 
In comparison to the otherwise more familiar random 
graphs studied by Erdos and Reny ||, |j, Bethe lattices 
at a given n and k avoid fluctuations in the connectivities 
of vertices and in the total number of bonds. 

There are slight variations in the generation of Bethe 
lattices. For instance, to add a bond one could choose 
at random two vertices of connectivities < k + 1 to link 
until all vertices are k + 1-connected. Instead, we have 
used the method described in Ref. (|] to generate these 
graphs. Here, all the terminals on the vertices form a 
list of n(k + 1) independent variables. For each added 
bond two available terminals are chosen at random to be 
linked and removed from the list. Furthermore, for al- 
gorithmic convenience, we reject graphs with possess self 
loops, bonds that connect two terminals of the same ver- 
tex. Multiple bonds between any pair of vertices are al- 
lowed, otherwise it is too hard to generate feasible graphs 
for, say, n = 32 and k+ 1 = 20. Since k + 1 remains finite 
for n —>■ oo, the energy and entropy per spin would only 
be effected to 0(l/n) by the differences between these 
choices. 

Once a graphical instance is generated, we assign ran- 
domly chosen but fixed couplings Jij £ {— 1, +1} to ex- 
isting bonds between neighboring vertices i and j. Each 
vertex i is occupied by a spin variable xi £ { — 1,+1}. 
The energy of the system is defined as the difference in 
number between violated bonds and satisfied bonds, 



H 



{borids} 



(i) 



and in this paper we will focus on the energy per spin 

H 



efc+i(n) 



(2) 



as a function of k + 1 in the limit of n — > oo. Each 
instance can have a large degeneracy VL in the configu- 
rations exhibiting its ground state energy, and we also 
sample the average entropy, 



s k+1 (n) = - lnO, 
n 



(3) 



III. r-EO ALGORITHM FOR BETHE LATTICES 

The extremal optimization algorithm, r-EO, which we 
employ in this paper, has been discussed previously in 
pl[ , an d m 0j |35| with regard to the setting of its one 
free parameter, r. Here, we merely describe the imple- 
mentation of r-EO without further justification. 

To obtain the numerical results in Sees. tV-^, we used 
the following implementation of r-EO: For a given spin 
configuration on a graph, assign to each spin Xi a "fit- 
ness" 



A, 



-^violated bonds = — 0, —1, —2, 



so that 



-(fc + l),(4) 



(5) 



is satisfied. Each spin falls into one of only k + 2 possible 
states. Say, currently there are spins with the worst 
fitness, A = — (fc + 1), rife with A = — k, and so on up to 
no spins with the best fitness A = 0. Now draw a "rank" 
I according to the distribution 



P(l) 



1 



1 - n 1 



:l~ r (l<l<n). 



(6) 



Then, determine < j < (k + 1) such that Yli=j+i n i < 



fc+i 



Finally, select any one of the nj spins in 



for these instances. 



state j and reverse its orientation unconditionally. As a 
result, it and its neighboring spins change their fitness. 
After all the effected A's and n's are reevaluated, a new 
spin is chosen for an update. 

This EO implementation updates spins with a (r- 
dependent) bias against poorly adapted spins on behalf 
of Eq. (^) . This process is "extremal" in the sense that 
it focuses on atypical variables, and it forms the basis 
of the EO method. The only adjustable parameter in 
this algorithm is the power-law exponent r. For r = 0, 
randomly selected spins get forced to update, resulting 
in a random walk through the configuration space which 
would yield poor results. For r — > oo, only spins in the 
worst state get updated which quickly traps the update 
process to a small region of the configuration space which 
may be far from a near-optimal solution. The arguments 
given in J35| and a few experiments indicate that r = 1.3 
is a good choice to find ground states efficiently on Bethe 
lattices. 

The algorithm never converges or "freezes" into a par- 
ticular state but perpetually explores new near-optimal 
configurations. It is, of course, easy to simply store the 
lowest energy state found so far in a given run of r-EO 
and terminate when desired. Previous experience with 
optimizing spin glasses with EO |2l| , and a few experi- 
ments, suggest a typical number of updates of 0(n 3 ) for 
an EO-run to obtain saturation in the values found for 
ground states, at least up to the system sizes n « 10 3 ob- 
tainable here. Instead of pushing to attain larger values 
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of n, we opt here for obtaining better statistics by sam- 
pling more instances at smaller values of n while spend- 
ing even more time on each instance than may seem 
to be required, in an attempt to ensure accuracy. In 
particular, our implementation restarts for each instance 



at least 



= 4 times with new random initial spin- 



assignments, executing « 0.1 n 3 updates per run. If a 
new, lower-than-previous energy state is encountered in 
run r, we adjust r max = 2 + 2r for that instance so that 
EO runs at least twice as many restarts as were necessary 
to find the lowest state in the first place. Especially for 
small n, r max hardly ever exceeds 4; for larger n a few 
graphs require up to 25 restarts before termination. 

Since EO perpetually explores new configurations it is 
well suited to explore also the degeneracy of low-energy 
states. In this case we not only store the first config- 
uration found with the lowest energy for that instance. 
Instead, we consider each configuration with the lowest 
energy, retaining new ones and rejecting all others. This 
procedure is somewhat inefficient and at best allows sys- 
tem sizes up to n — 256 beyond which the degeneracy 
exceeds memory constraints. But it provides a fast way 
to also determine the T = entropy of the ground states 
with moderate accuracy. In these runs, we used a similar 
approach to the above, except for setting r max = 8 + 2r 
where r is the latest run in which another new configu- 
ration of the lowest energy was located. Here, for some 
highly degenerate instances at larger n, r max could reach 
up into the 100's, further limiting attainable system sizes. 



IV. NUMERICAL TEST 

To evaluate the proposed r-EO algorithm, we have run 
a series of test. First, we can defer to some already pub- 
lished results |]2lL p5| . In Ref. Q we have calculate 
approximations to the ground state energy for ± J spin 
glasses on a hypercubic lattice for d — 3 and d = 4 for sys- 
tems up to n = 12 3 = 1728 which for each n reproduced 
previous results obtained with sophisticated genetic al- 
gorithms |27| ^8| (although there we used a fixed r max ). 
To evaluate the ability of the algorithm to determine the 
degeneracy of low-energy states found, we have repro- 
duced within statistical error the results of Ref. |36 up 
to n — 6 3 beyond which EO ran out of time and memory 
to sample states completely. (Ref. |36| used a more effi- 



cient way to estimate the entropy from sampling only a 
small number of states.) And it took EO only a fraction 
of a second to find all 60 ground states of a 4 3 instance 
that had been exactly enumerated in Ref. f23fl . 

To gauge r-EO's performance for larger n, we have 
run our implementation also on two 3d lattice instances, 
toruspm3-8-50 and toruspm3-15-50, with n — 8 3 = 512 
and n = 15 3 = 3375, considered in the 7th DIMACS 
challenge for semi-definite problems |?7| . Bounds on 
the ground-state cost established for the larger instance 
are H\ OWCI = —6138.02 (from semi-definite programming) 
and -ffuppcr = —5831 (from branch-and-cut). EO found 



H EO = -6049 (or H/n = -1.7923), a significant im- 
provement on the upper bound and already lower than 
lim^oo H/n w 1.786 . . . found in Refs. JHJ [27], |§. Fur- 
thermore, we collected 10 5 such states, which roughly 
segregate into 3 clusters with a mutual Hamming dis- 
tance of at least 100 distinct spins; at best a small sample 
of the w 10 7 3 ground states expected [^6j! For the smaller 
instance the bounds given are -922 and -912, resp., while 
EO finds -916 (or H/n = —1.7891) and was terminated 
after finding 10 5 such states. While this run (including 
sampling degenerate states) took only a few minutes of 
CPU (at 800MHz), the results for the larger instance re- 
quired about 16 hours. 

Finally, we note that we have considered the algorithm 
for making Bcthc lattices previously in Refs. 31, |35[ . 
In Ref. Iplfl we have studied the graph bipartitioning 
problem and found that the ground state energy was 
well above previous RS calculations from Ref. [§|, but 
only minutely below numerical calculations obtained us- 
ing simulated annealing jl9| . In Ref. [[35] we have consid- 
ered some variations in the generation of Bethe lattices 
and found that they effect the results only in next-to- 
lcading order. 



V. NUMERICAL RESULTS FOR BETHE 
LATTICES 

We have simulated Bethe lattices with the algorithm 
described in Sec. H for k + 1 between 3 and 26, and 
graph sizes n = Irfor I = 5,6,... ,10 to obtain re- 
sults for ground state energies, and for n G [16 . . . 256] 
to determine their entropy. In the following, we present 
the results for ground-state energies and entropies from 
those simulations. The results are discussed in detail in 



Sec. VC 



A. Ground State Energies 

To reach relative statistical errors of our averages 
roughly uniform with n we generated initially a number 
of 10 5 /y/n instances for each n and k + 1. Fortunately, 
deviations appear to narrow much faster than 1/y/n, and 
thus we added more instances at smaller n with small ex- 
tra computational cost to obtain narrow error bars there 
as well. In Tab. Ejwe list the values of average energies 
according to Eq. (g), (e fc+1 (n)), for each k + 1 and n, the 
number of instances used, and the average number of up- 
date steps required. Tab. [n] lists a few properties of the 
computations. The results for the number of updates has 
been also averaged over all connectivities + although 
lower-connected graphs require typically fewer updates. 
Note that this is the minimal number of updates needed 
to obtain the listed results, the actual number of updates 
taken up by each run of EO to ensure convergence was at 
least twice of that but could be much larger, according 
to the specification of the algorithm in Sec. III. 
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TABLE I: Data from the EO simulations for the average ground-state energy per spin ek+i(n), plotted also in Figs. |l]and ^. 



n — eg(n) — ej(n) — eg(n) — gg(n) — g7(n) — es(n) 



-eg(n) 



-eio(n) — ei 5 (n) — e 2 o(n) — e 2 s(n) 



32 1.3506(6 
64 1.2231(3 
128 1.2426(5 
256 1.2542(3 
512 1.2608(2 
1024 1.2644(1 
2048 1.2673(1 



1.5543(6) 1.6734(8) 1.8424(8) 1.9425(9) 2.0906(9) 2.1730(33) 2.7013(40) 3.1056(48) 3 

1.3964(4) 1.5979(4) 1.7294(5) 1.8972(5) 2.0083(6) 2.1551(6) 2.2557(23) 2.7884(28) 3 

1.4245(10) 1.6269(10) 1.7652(12) 1.9335(13) 2.0476(14) 2.1945(15) 2.3042(16) 2.8443(20) 3 

1.4417(6) 1.6434(7) 1.7885(8) 1.9549(8) 2.0782(10) 2.2204(10) 2.3324(11) 2.8774(14) 3 

1.4534(4) 1.6548(5) 1.8020(5) 1.9685(6) 2.0934(6) 2.2379(7) 2.3488(7) 2.8993(10) 3 

1.4603(3) 1.6612(3) 1.8110(3) 1.9762(5) 2.1035(5) 2.2470(5) 2.3605(5) 2.9092(7) 3 



4923(51) 

2090(33) 3.6029(37) 

2742(23) 3.6731(25) 

3186(16) 3.7215(18) 

3435(11) 3.7505(12) 

3551(9) 3.7612(11) 



oo 1.2719(5 



1.472(1) 1.673(1) 1.826(1) 1.990(3) 2.121(1) 2.2645(5) 2.378(3) 2.935(1) 3.389(1) 3.806(4) 



TABLE II: Some properties of the numerical computations. 
Listed are for each n the number of instance used and the 
average number of updates for each instance needed to obtain 
the results listed in Tab. M. 



n 


Instances 


t 


32 


19444 


3.0 10 2 


64 


13750 


1.5 10 3 


128 


883 


1.0 10 4 


256 


625 


1.610 s 


512 


441 


3.1 10 6 


1024 


312 


7.6 10 7 


2048 


220 


1.5 10 8 



Unfortunately, when plotted as a function of 1/n, the 
average energies for each given k+1 clearly do not extrap- 
olate linearly (as, for example, seems to be the case for 
spin glasses on a hypercubic lattice [pl"[ 2?] , p8|). Instead, 
using an extrapolation according to K5 1 



e fe+ i(n) ~ e fe+ i + ^- (n — > oo). 



(7) 



We find that for the whole range of connectivities k + 
1 studied here, the scaling corrections appeared to be 
consistent with v — 2/3 within a few percent, except 
for two outliers at k + 1 = 10 and 25. Thus, we have 
plotted for each k+1 the values of ek+i (n) as a function 
of 1/n 2 / 3 in Figs. [I] and [| Although the extrapolation 
appears to be linear on that scale for each k + 1, we have 
fitted the data with the more general form of Eq. (0). 
(Fits were weighted according to n and to the inverse of 
the standard deviation for each point.) These fits are 
also shown as dashed lines in each of the Figs. |l| and ||. 

The extrapolation results for the ground state energies 
appear to be quite stable under variation of the scaling 
form, for instance, when fitting with fixed v = 2/3 in- 
stead of Eq. (Q). We estimate that each has a relative 
error of about jO.3%. Exceptions to this estimate we 
have to grant for the cases of k + 1 = 10 and 25, in which 
case we also observe significant differences to the v = 2/3 
corrections to scaling. 

We can compare our results with existing theoretical 
predictions at the RS and the 1RSB level at least for the 




FIG. 1: Extrapolation plot for the EO data for k+1 = 3 given 
in Tab. | and the fitted curve according to Eq. (Q). (Data 
points were weighted with respect to n and the inverse of the 
error.) For n — > oo the extrapolation gives ez = —1.2719(5), 
way above the RS result but consistent with the 1RSB result 
from Ref. |hj, both indicated by horizontal lines. 



case of k + 1 = 3. For this case, a recently published 
calculation yielded — —1.2777 at the RS, and e-$ = 
— 1.2717 at the 1RSB level (further replica corrections are 
expected to be small). These values are also indicated in 
Fig. 0. Clearly, our result for k + 1 is consistent with 
the 1RSB results, but certainly inconsistent with the RS 
result. Further 1RSB results for other values k + 1 are 
currently being calculated (3!|. We will discuss a more 
detailed analysis of the extrapolated values of ek+i at 
n — > oo in Sec. VC. 



B. Ground State Entropy 

We have also used EO to sample the degeneracy 51 
of the lowest-energy states found. Due to the discrete 
nature of the energy of the system, ground-states can be 
highly degenerate, and the ground-state entropy per spin 
defined in Eq. (^J) may well be non- vanishing for n — > oo. 
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FIG. 2: Extrapolation plot for the EO data in Tab. [j] for k + 1 = 4 to 25, as in Fie. |IJ All data seems to extrapolate well 
linearly in l/n 2 ^ 3 . The extrapolated values of e^+i for n — > oo are also listed in Tab. 
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While the search for a ground state of an instance is 
certain to provide a rigorous upper bound to the actual 
ground state energy, the search for the complete set of 
ground states for an instance entails the risk of a two 
competing systematic errors. (1) If EO misses to find 
the exact ground state, one is likely to vastly over-count 
the degeneracy, since is expected to rise exponentially 
with the energy above the ground state po| . (2) Even 
if EO finds ground states, it may simply undercount O, 
since such states could be too far separated in configura- 
tion space. Therefore, we have implemented EO with the 
settings described in Sec. Ill , which emphasize the desire 
for accuracy over computational efficiency. Accordingly, 
we were bound to conduct a separate set of simulations 
from those that determined the energies only. In these 
simulations we focused on smaller system sizes of n < 256 
for k + 1 = 3, . . . , 9 and 10, 14, . . . , 26 only. The limit 
on n for the smaller k + 1 is mostly dictated by avoiding 
system sizes at which SI typically exceeds 10 6 . 

As a test for the accuracy of our implementation, we 
have run the simulation for k + 1 = 3 twice on the exactly 
identical instances, using different initial conditions and 
n/5 more updates in the second run: The results, both 
for the energies and f2, were identical for each instance, 
producing the same set of configurations independent of 
the starting point of the search. We therefore assume 
that systematic errors in our data are small and can be 
neglected. 

Since the range of system sizes n is smaller than for the 
case of the energies, it is more difficult to extrapolate our 
data for (sk+i(n)). Again, it is clear that the corrections 
are not linear in 1/n, but instead seems to be scaling 
close to 1/n 2 / 3 for all k + 1, as for the energies above. 
Considering the limitations on n, we assume that the 
corrections are exactly of that form and extrapolate our 
data simply with a fit to 



Sk+i(n) ~ s k+ i + 



A 

7m 



(n — > oo), 



(8) 



again, weighting each data point with respect to n and 
the inverse of its error. While the systematic and sta- 
tistical uncertainties of our data appears to be small, 
the uncertainty about the scaling corrections must be 
considered the most significant limitation on accuracy in 
our extrapolation. The data and the extrapolation fits 
according to Eq. (||) are shown in Figs. || and [l| The 
results for s^+i for n — > oo are listed in Tab. HJ. 

The data clearly shows a different quantitative behav- 
ior between odd and even values of k + 1. This differ- 
ence for the entropies can be explained in terms of the 
"free spins:" In a highly frustrated system, even near 
ground states, many spins are stuck in a situation in 
which they violate many of their constraints, no mat- 
ter how they are oriented, and changing from one direc- 
tion to the other may hardly change the energy of the 
system. In particular, an even-connected spin that hap- 
pens to violate exactly half of its bonds (with J = ±1) 
can flip freely without any change to the energy. Odd- 



TABLE III: Extrapolation results for the entropies per spin 
for the data plotted in Figs. |3 and [l|. 



k+ 1 


Sfe+l 


k + 1 


Sfc+l 


3 


0.0102(10) 


4 


0.0381(15) 


5 


0.0048(10) 


6 


0.0291(10) 


7 


0.0020(10) 


8 


0.0218(10) 


9 


0.0002(15) 


10 


0.0198(10) 


15 


0.0002(15) 


14 


0.0126(10) 






18 


0.0095(10) 






22 


0.0076(10) 






26 


0.0063(15) 



connected spins can only become "free" in a connected 
pair (that happens to violate exactly half of its exter- 
nal bonds but satisfies their mutual bond) in which both 
simultaneously flip without changing the energy. The lat- 
ter situation is naturally far less likely, and thus, purely 
even-connected graphs exhibit far more potential for de- 
generacy at the ground state than the corresponding odd- 
connected graphs. Some preliminary studies for k + 1 = 3 
and 4 show that in ground state configurations the frac- 
tion of free spins (zero by design for k + 1 =3) converges 
to a value just around 5% for k+ 1 = 4, while the fraction 
of free pairs seems to vanish for large n for both, even 
and odd k + 1. We have not explored the clustering of 
these states |$6|. We will explore the different behaviors 
for even and odd k + 1 in the next Section. 



C. Discussion of the Extrapolation Results 

In this section, we want to focus on some of the curi- 
ous properties exhibited by the values of the energies and 
entropies found by extrapolation in the previous section. 
We have already noted the difference between the en- 
tropies for even and odd values of k+ 1. In fact, there are 
similar differences, although more subtle, for the energies 
ek+i- These differences become most apparent when we 
plot the data asymptotically for large k + 1, where it is 
known that 



lim 



efc+i 



k+l^oo Vfc + 1 



E SK , 



(9) 



with Esk — 0.7633 being the RSB ground state energy 
of the Sherrington-Kirkpatrick model || [ll). In Fig. || 
we have plotted ek+i/yk + l as a function of l/(k + 1). 
On this scale, we notice that the energies split into a set 
of even and a set of odd values, each located apparently 
on a straight line. Even though k + 1 < 25 is quite small, 
each line separately extrapolates very close to the exact 
value for large k + 1 indeed: Egg » -0.763 and Egg m 
—0.765. Even more amazing, the value of ei = —1 [see 
Eq. ( |l6| ) below] for the trivial k + 1 = 2 Bethe lattice 
is very close to the linear fit for the even EO results. 
Clearly, a function that would interpolate continuously 
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FIG. 3: Extrapolation plot for the EO data for the entropy Sk+i(n) for k + 1 = 3 to 15. All data seems to extrapolate well 
linearly in l/n 2//3 . Note the difference in the results between odd (left) and even (right) k + 1. The extrapolated values of Sk+i 
for n — * oo are listed in Tab. III. 
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FIG. 4: Extrapolation plot for the EO data for the entropy 
Sfc + i(n) for some larger, even k + 1, similar to Figs. H. 



all the data will have to be very complicated (oscillatory). 
But could it be that its envelope on the even and the odd 
integers happens to be simple? Then, in case of the even 
data [H , we could even write down the exact form of the 
function for £k+i that would fit the data, since we know 
it also has to pass e2 = — 1 and satisfy Eq. 



£k+i = Vk+lEsK 



2E SK + V2 



s/k + l 



(10) 



To test Eq. (|1C|), we plot the data in Fig. |6| as 
ek+\/ £k+i to study its deviations from the conjecture. 
While the extrapolated values do not fall exactly within 
their (estimated) error bars on the proposed form, they 
are indeed within about 0.1% of it. To judge how close 
the data is to the proposed functional form in Eq. (10), 
we utilize a closely related example. The ground-state 
energy as a function of the (continuous) average connec- 
tivity (c) is known exactly for the RS case of ordinary 
random graphs with fluctuating internal connectivities, 
Eq. (16) in Ref. p2[. If one plots that solution (which 
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FIG. 5: Plot of the rescaled extrapolated energies, 
Ck+i/Vk + 1, as a function of l/(k + 1). The data appears 
to fall on two separate straight lines for even and for odd 
k + 1. The straight line provides an excellent fit all the way 
from the exact result e2 = — 1 (diamond) to Esk = —0.7633 
(horizontal line) at k + 1 — > oo. 



involves exponentials and modified Bessel functions) in 
the same way as ek+i in Fig. one notes that it, too, 
could be approximated surprisingly well with a straight 

line, —yj2c/n + (■y t /2/n — |J /s/c, now crossing the RS 

ground state energy —y/2pK for (c) — ► oo and reach- 
ing the trivial result of —1/2 at the percolation point 
(c) = 1. In Fig. ^ we superimpose the relative error of 
this approximation with respect to the exact RS result 
with the relative error of our data with respect to the 
conjecture. It shows that the error of the conjecture is 
still almost by an order of magnitude smaller than the 
global bound for the RS example, thus putting a signifi- 
cant bound on any corrections similar in type to Eq. (16) 
in Ref. ||. 

The differences between even and odd connectivities 
are even more pronounced in case of the entropies, as 
we have explained in Sec. [v|. Thus, although our data 
for the entropy is not nearly as accurate as for the ener- 
gies, it is still instructive to study it in more detail. In 
Fig. ||, we plot the extrapolated values of the entropies 
from Tab. Ill to explore its decrease for large k + 1. De- 
spite the large error bars, a significant qualitative differ- 
ence between even and odd data points is visible: The 
entropy for even values of k + 1 decays slowly, apparently 
linearly with l/(fc+ 1). On the other hand, the entropies 
for odd fc + 1 drop much more rapidly, and are already in- 
distinguishable from zero (within our errors) for fc+1 = 9, 
while it is clearly non- vanishing for k + 1 = 3 [unless our 
assumption about the scaling corrections in Eq. (H) are 
incredibly wrong (see Fig. ^)]. Unfortunately, with only 
a small, discrete number of data points available that are 
significantly above zero, it is very hard to decide whether 
the entropy for odd fc+1 merely decays exponentially, or 
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FIG. 6: Plot of the energies relative to the conjectured func- 
tion £k+i in Eq. ( |lo| ) as a function of l/(k + 1). All data for 
even k + 1 falls within about 0.1% of Ek+i (i- e. the horizontal 
line). The point at k + 1 — 2 (diamond) is exact by definition, 
of course. 



FIG. 8: Asymptotic plot of the extrapolated entropies from 
Tab. £n| as a function of l/(k + 1). The data for even k + 1 
seems to vanish linearly with 1 / (k+1) (dashed line) . The data 
for odd k+1 drops more precipitously, and can not reasonably 
be fitted at this level of accuracy. 
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FIG. 7: Plot of the relative error e in % between the extra 
olated data and the function in Eq. (|lC|), as before in Fig, 
but now superimposed with the corresponding result obtained 
for the exactly solvable RS spin glass on random graphs [ {lif . 
Note the oder-of-magnitude larger deviations from the refer- 
ence line for the RS example. 



whether there exists a finite value of k + 1 above which 
all odd entropies become identically zero. 



VI. CONCLUSION 

In this paper we have presented a extensive numerical 
study of the ground states of spin glasses on Bethe lat- 
tices. The available data possessed sufficient accuracy to 
obtain extrapolated values for ground state energies and 
entropies at the 0.1% and the 10% level, respectively. In 



both cases, we significant differences 
the data for odd and even values of 
the numerical results, we showed that 
energies for all even values of 2 < k -\ 
fitted with a simple function, Eq. (|io|). 
data suggests that the entropies are g 
at small k + 1, but may vanish above 
odd values. 



emerged between 
k+1. Based on 
the extrapolated 
- 1 < oo was well 
Furthermore, the 
enerally non-zero 
a finite k + 1 for 



Of course, there is plenty of reason to doubt that such 
a simple result as Eq. (|I(]), albeit confined to discrete 
integer values of k + 1, could indeed be the solution to 
a complex RSB problem. In fact, one argument against 
the conjecture is a discrepancy in its prediction for large 
k + 1 at next-to-leading order. Several authors [||, [T^] 
have studied spin glasses on random graphs beyond the 
RS level perturbatively for k + 1 = z — > oo to determine 
the 1/z correction /i to the free energy (at T > 0) in 
Esk + fi/z. It was recently calculated ||, that the cor- 
rection for fixed connectivities for T — > would be about 
/i = -0.317, while Eq. © would predict about 0.1124. 
It should be noted, though, that the 1/z expansion im- 
plicitly assumes a smooth continuation off the integers 
which may lead to ambiguities in light of the oscillatory 
behavior between even and odd integers we found for the 
ground states in Fig. || (similar to the continuation of, 
say, the function cos(27rz)/z for z — ► oo, although this 
function would not possess a 1/z expansion at all). 

In any case, future calculations like the one in Ref. 
but for even k + 1 will provide a check on both, our 
extrapolated data and the conjecture. 
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where i > 0. It is easy to show that these equations 
satisfy the constraints in Eqs. (fi"l|). 

We can transform these equations by defining 9 = 
t/n, dd = 1/n, y(x, 9) = ^Y^*L k,tx\ and p(x,6) = 



i=0 I 



Considering n large and 9 continuous, 



Eqs. (|12| ) turn into 



Appendix 

A. The Case k + 1 = 2 

Clearly, a Bethe lattice in which each vertex has ex- 
actly 2 connections can only consists of a collection of 
disconnected loop graphs. We merely need to determine 
the number of loops and their size distribution to derive 
the average ground state energy and entropy. Each loop 
has a 50% chance of being frustrated, thus, the number 
of cut bonds is equal to one-half of the number of loops, 
and the degeneracy is equal to the length of these loops 
to the power of their number. 

To analyze the k + 1 = 2 case we consider each vertex 
as a node with two terminals. Adding lines can create 
two types of objects: strings and loops. We consider, 
after adding t lines, an individual vertex as a string of 
length 0, of which there are Zo,t; in general, we have 
strings of length i, each possessing two open terminals. 
In particular, before we added any lines: h.t=o = $i,o- A 
loop of length i is created by addition of a line to both 
open terminals of a string of length i — 1. There are 
Pi t loops of length i after adding t lines which can not 
evolve further, since they don't possess any more open 
terminals. We start with 2n open terminals and cover 2 
of those with each newly added line. We can identify two 
constraints: 



oo 

E 

i=0 



E 

i=l 



i (h,t +Pi,t) = t. 



(11) 



After adding t lines at random, there are 2(n — t) ter- 
minals left to accommodate the next line, allowing for 
(2{n-t)-^ (jjfj eren ^ choices. Accounting for all possible 
choices, we obtain 



h,t+l 



Pi,t+1 = 




dy(x,6) 
d9 



2y(x,6) 



x[y(x,0)? 



1 - 



(l-0)[l-0-l/(2n)] 



+- 



n 



y(x,9) - xy(x 2 ,9) 



2(l-0)[l-0-l/(2n)] 



dp(x,9) xy(x,9) 
d9 = (l-6>) 2 ' 

y(x,0) = l, p(x,0)=0. (13) 

Luckily, for n — > oo, the equations are easily solved to 
give 



p(x,9) = --\n(l-x9). (14) 



Finally, the total number of loops for the (almost) com- 
pleted graph, 9 = 1 — 1/n, is given by 



oo 



1 



p(l, 1 - 1/n) = Pi >n -l ~ g ln ( n )' 
»=l 



(15) 



On average, half of these loops will be frustrated, i. e., 
they will have one of their bonds violated. Since the 
Hamiltonian in Eq. ( ^) counts the difference between vi- 
olated and satisfied bonds, or twice the violated bonds 
minus the number of all bonds, n(k + l)/2 = n, we get 



H ln(n) 

e= i + ^r^- 

n 2n 



(16) 



Similarly, we can calculate the degeneracy f2 of these 
ground states, roughly, as the average length of loops, 
(i) = d x lap(x, 1 — l/n)|a; = i ~ n/ln(n), taken to the 
power of one-half of their number, ln(n)/(4n), to give 



s = — In il 

n 



ln(n) 2 
4n 



(17) 



Clearly, both the number of violated bonds as well as the 
entropy vanish in the large-n limit. 
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